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ABSTRACT 



Solar activity and helioseismology show the limitation of the standard solar model and call 
for the inclusion of dynamical processes in both convective and radiative zones. In this paper, 
we concentrate on the radiative zone. We first recall the sensitivity of boron neutrinos to the 
microscopic physics included in solar standard and seismic models. We confront the neutrino 
predictions of the seismic model to all the detected neutrino fluxes. Then we compute new 
models of the Sun including a detailed transport of angular momentum and chemicals due to 
internal rotation that includes meridional circulation and shear induced turbulence. We use two 
stellar evolution codes: CESAM and STAREVOL to estimate the different terms. We follow 
three temporal evolutions of the internal rotation which differ by their initial conditions: very 
slow, moderate and fast rotation, with magnetic braking at the arrival on the main sequence 
for the last two. bf We find that the meridional velocities in the present solar radiative zone 
is extremely small in comparison with those of the convective zone (smaller than 10~^ cm/s 
instead of m/s). All models lead to a radial differential rotation profile in the radiative zone but 
with a significantly different contrast. We compare these profiles to the presumed solar internal 
rotation and show that if meridional circulation and shear turbulence were the only mechanisms 
transporting angular momentum within the Sun, a rather slow rotation in the young Sun is fa- 
vored. We confirm the small influence of the transport by rotation on the sound speed profile but 
its potential impact on the chemicals in the transition region between radiation and convective 
zones. These models are physically more representative of the real Sun than the standard or 
seismic solar models but a high initial rotation, as it has been considered previously, increases 
the disagreement with neutrinos and the sound speed in the radiative zone. This present work 
pushes us to pursue the inclusion of the other dynamical processes to better reproduce the ob- 
served solar profile in the whole radiative zone and to better describe the young active Sun. We 
also need to get a better knowledge of solar gravity mode splittings to use their constraints. 
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Subject headings: hydrodynamics - neutrinos- Sun : evolution - Sun: helioseismology -Sun 
: interior -Sun: rotation 



1. INTRODUCTION 



Evidence for the presence of dynamical processes exist in stars all over the Hertzsprung-Russell dia- 
gram. In the case of the Sun, rotation and magnetic activity have been confirmed for more than 4 centuries. 
Rapid rotation is found mainly in young stars and in termediate-m ass to massive stars while mild to slow 
rotation is found in low-mass stars and in giant stars. iKrafti (|l965h was the first to study in details the pro- 
jected rotational velocity vsinj of young stars in clusters and shows the transition between rapid rotation for 
early-type stars and slow rotation for late-type stars. Then, Weber explained that the angular mo mentum loss 



is due to a magnetized wind in solar-like stars. In studying Pleaides, Ursa Major and Hyades, ISkumanich 



(Il972h used solar and cluster data to establish the law of variation of the rotatio n with age, this simple law is 
in fact more complex and has been studied in great d etails for different masses Jstauffer et al.l [l98 71: IB ou vied 



19941 : iMaeder & Mevnetll200d 120041 : lBouviedl2009h . 



Rotation affects the internal stellar structure and evolution both directly via the modification of the 
gravitational potential, and by means of associated transport processes. It has a direct impact on the shape 
of stai's, that can be directly probed with interferometric data for the sti'ong l y deformed fast rotating stars 



(Domiciano de Souza et al. 



2006 



2OO3I: iMcAlister et alboosl : Ivan Belle et albood : 



Kervella & Domiciano de Sousa 



Domiciano de Souzal2008r) . In the case of the Sun, the effect is small as shown in 



Piau & Turck-Chieze 



(l2002h but of considerable importance (Zahn 2009; Rozelot 200 9) to check the influence of the deep core 
rotation and of the deep magnetic field (IDuez et al.ll2009al . l2010l) on the tachocline and on the solar surface. 



The standard solai^ and stellar evolution models do not take into account the effects of such dynamical 
processes as rotation or magnetic fields. If helioseismology required the improvement of the solar model 
and motivated the introduction of microscopic diffusion, going beyond the zero-order model, it is the far 
too rough agreement between observational data and theoretical predictions for other stars that has led to 
progressively inti'oduce the dynamical physical processes likely to affect the trans port of momentum an d 
chemicals in the stellar evolution codes. Using e i ther a purely d iffusive approach (Endal & Sofial (119811) : 
Pinsonneault et al. I (ll989l):IChabover et al.l(ll995h:lLangeii (ll998l) : lHeger et all (|2000l) ) or the more complex 
formalism developed by ZahnI ( 19921 ) and Maeder & Zahn (1998), the introduction of the rotation-induced 
transport of angular momentum and c hemicals to model both massive and low-mass star s improves the 



comparison with observations (see e.g. iMeynet & Maedenl2000l: iTalon & Charbonnellll998l : IPalacios et al 



2003h . These results encourage us to pursue this effort for the Sun to get a proper interpretation of the 
existing helioseismic observations and the coming asteroseismic ones. 



'P.A.P Nghiem has now left the Service d'Astrophysique in IRFU/CEA and works in another discipline in the same institut. 
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Acoustic and gravity mo de detections provide a unique insight on the inte r nal solar sound speed , 



density and rotation profiles (IKosovichev et al 



Couvidat et alJl2003bl : iTurck-Chieze et alJl2004al : 



1997 



Mathur et al 



2007 



Turck-Chieze et alj 12001 al : [Thompson et all l2003l : 



20081) . This gives us a unique oppor- 



tunity to validate the complex formalism introduced in rotating stellar evolution models. This information 
justifies the development of a Dynamical Solar Model to be confronted with the helioseismic and neutrino 
probes. The solar status must improve our knowledge of solai^-like stars where only external stellai^ rotation 
rates or abundance anomalies are used to confront theoretical assumptions to observational facts. Moreover, 
building a consistent dynamical evolution model for the Sun will largely contribute to establish a complete 
and consistent MHD representation of the Sun and a good connection between internal and external mag- 
netism. 

The inner solar rotation is the first evidence of the internal dynamics that needs to be understood. 
This is not an easy task because the present-day profile results from the interplay between several dis tinct 
processes dZahn et al.ll 19971 : IChaboyer et al.lll995uEggenberger et al.ll2005MCharbonnel & Talonll2005l) and 
their influence depends on the adopted theoretical prescriptions. The most recent of these studies include the 
transport of angular momentum by magnetic field or internal gravity waves, in addition to the "purely" rota- 
tional transport by meridional circulation and shear-induced turbulence, to obtain a more efficient extraction 
of angular- momentum and a relative flat rotation profile in the radiative interior in better agreement with 
helioseismic data. Actually, both magnetic field and internal gravity waves probably contribute to the trans- 
port of angular momentum in stellar interiors, yet no models have been published up to now including all 
the processes. Moreover, the way to account for the effect of magnetic fields in 1-D stellar evolution codes 
is still puzzling and a matter of debate, and the introduction of grav ity waves produces some irregularities 
in the radiative rotation profile (see e.g. Icharbonnel & Talonlboosi) in disagreement with the heUoseismic 
results. 



In the present paper, we limit ourselves to the sole effect of rotational transport and carefully study the 
hypotheses and the order of magnit ude of the terms that we introduce for the inner ro t ation. This approach 
pushe s further the previous works bvlChabover & ZahnI (119921) : IChaboyer et al.l(ll995r ): ICharbonnel & Talon 
(|2005l ): lEggenberger et al.l (120050 : 1 Yang & Bil(|2006l.l2007l ). We focus our analysis on the solar core and we 
compare our results with all the recent existing seismic and neutrino indicators. This paper is a first of a 
series where we will discuss in details the influence of rotation, magnetic field and internal waves on the 
solar dynamical model. We do not present and disc uss the abundances of lithium and beryllium as it has 
been done extensivelv bv IPinsonneault et al. I (|l989l) . Indeed we have shown in our previous works on the 
tachocline (BrunetaL 19991) and on young stars that the abundance of these elements are sensitive to the 
early evolution as well as to the presence of the tachocline. These studies lead us to the conclusion that we 
need to tre at properly the magnetic fie ld in the early phase to better estimate the lithium evolution during 
this phase (IPiau & Turck-Chiezell2002l) . 



In Section 2, we present the status of the solar classical models and compare the seismic mode l pre- 



Zahn 



diction s to all the present de t ected neutrinos. Then in Section 3 we recall the formali sm proposed by 
(|l992l ) and Maeder & ZahnI (|l998r) . and shghtly improved by Mathis & ZahnI \20o4) that we use to intro- 
duce the rotation effects in the stellar evolution codes. In that section, we also give a brief description of 
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the CESAM and STAREVOL stellar evolution codes used in the paper, and present the three different types 
of solar models which differ by their initial rotation and the presence of magnetic braking. In Section 4 
we compare the results of these three rotating solar models obtained with both codes. These models are 
compared with seismic observations and neutrino detections in Section 5. In the last section, we summarize 
the important points in a more general context. 



2. THE MICROSCOPIC PROCESSES AND THE SOLAR PROBES 

The Sun differs from other stars by its proximity and our ability of observation. We know better its 
luminosity, radius and mass than for any other star. This has, for long, turned the Sun into a reference for 
stellar evolution. Nevertheless the numerous constraints point to flagrant limitations of the standard picture: 



In this framework, the luminosity increases by 30% in 4.6 Gyr, 8% during the last Gyr and by only 
10"^ during th e last century. B ut in fact, the so-called solar constant shows variations by about 10"^ 
in mean value (IFrohlichlbOOq) coiTclated to the solai^ cycles (ACRIM and SoHO observations). The 
presence of faculae around the sunspots explains the increase of luminosity at the maximum of the 
cycle. Moreover the series of data shows clearly short time scale luminosity variation connected 
to the rotation and the presence or absence of sunspots in the field of view. This fact shows, just 
from fundamental quantities, that the standard framework is not sufficient and that magnetism and 
rotation must be introduced in the solar structural equation s to interpret these observations because it is 



believed that these manifestations have an internal origin (jTurck-Chieze & Lambertll2007l : iDuez et al. 
2009ah . 



We have learned recently that the absolu te irradiance is no t so well established. The SORCE satellite 
suggests a reduction of nearly 6 W/m^ (|Kopp et al.ll2005|) in comparison with the standai^d (3.846 ± 
0.00410^^^ ergs/s corresponding to 1367.6 W/m^) used for solar model, see Bahcall & Pinsonneault 
(1995). This result modifies by 1.5-2 W/m^ (340 instead 342 W/m^) the energy which reaches the 
stratosphere. So, the present s olar luminosity value, which is u sed to calibrate the standard model, 
is presently uncertain at 0.5% (|Turck-Chieze & Lefebvrell2010h . The calibration and ageing of the 
radiometers are difficult to estimate, so ffie use of the space data is not easy. The SORCE/TIM in- 
strument has an estimated absolute accuracy of 350 ppm and a long-term repeatability of 10 ppm per 
year. Never theless, the ana lysis of the 30 years data series seems to show a slow decrease of ffie total 
luminosity (IFrohlichl 120061) between successive minima which cannot be understood in the standard 
solai- framework. 

For three decades, our star is also being scrutinized by two probes of the interior (neutrinos and 
seismic modes) that help us to largely progress on our capability to check the internal solar plasma 
in great details. We are now able to disentangle the production of neutrinos for specific nuclear 
reactions like ^Be(p,y) and ^B(p, y), and to sum the different flavours to really estimate the number 
of emitted neutrinos. The quality of the detection of the boron neutrinos (see § 2.2 and Tab. [B, 
is such that the central temperature of the Sun is now known to about 0.5%. Neutrinos can ffius 
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bring very strong constraints to the solar models. On the helioseismic side, the detection of million 
acoustic modes stimulated a real insight on the thermodynamics of the radiative region (see below) 
and on the dynamics of the convective zone. Gravity modes appear also promising to reveal the last 



missing information, the dynam ics of the solar core (ITurck-Chieze et al. 



2004bl : iGarcia et all 12007 . 



2008; Mathur et al. 2007, 2008). 



2.1. The validity of the standard solar model 



The solar radiative zone represents 98% of the total mass of the Sun and the adopted microphysics 
plays a basic role in this region. The equilibrium between gravitational energy, nuclear energy produc- 
tion and the energy escaping by photon interaction needs to be followed in the radiative interior over long 
time-scales. This region is now permanently probed by helioseismology, which is a key to validate the 
various ingredients used in the construction of the standard solar model. One success of the SoHO space 
observatory has been obtained by measuring the Doppler velocity shifts down to the low frequency range 
of the acoustic spec trum with two instruments: GOLF, Global Oscil lations at Low Freque ncy described by 
Gabriel et al.l (119951 ) and MDL a Michelson Doppler Imager built by IScherrer et al.l (|l995n . This part of the 
spectrum, contrary to the high frequency range, i s not sensitive to the v ariability of the sub surface layers 



along the solar cycle (ITurck-Chieze et al 



2001a 



Couvidat et al 



2003bl). The corresponding modes have 



longer lifetimes but smaller intensities (IBertello et al.ll2000l : lGarcfa et al.ll200lh . From these observations, 
a very clean sound speed profile has been established down to 0.06 Rq together with a reasonable density 
profile. It has been compared to the theoretical models for different solar compositions For t he latest sets 
of so lar abundances with lowered oxygen abundance (IHolweged l200ll : lAsplund et al.l |2009|; ICaffau et al. 
l2008k the theore tical profiles show clear di f ference s compared to the observed soundspeed, that still need 
to be understood ( Turck-Chieze et al.ll2001al.l2004a ). 



2.1.1. The nuclear processes 



It is interesting to notice that each physical ingredient of the structural equations (specific nuclear rate, 
specific opaci ty coefficient, screening o r Maxwellian tail distribution) has a specific signature on the sound 



speed profile (ITurck-Chieze 



et al 



1997h. It has t hus bee n of prime importance to check the validity of the 



involved nuclear processes. iTurck-Chieze et al.l (|2001br) note that the present sound speed profile does not 
favour any tiny variation of the Maxwellian distribution of the ions nor strong screening or large mixing 
in the core. The core sound speed profile and the specific signature of the pp reaction rate put strong 
observational constraint on the cross section at a level of 1% including the screening effect. This cross 
section is known only theoretically due to the weak character of the interaction. The other nuclear cross 
sections have been measured in laboratory during the last three decades but the extrapolation towards the 
stellar plasma conditions has been only measured for the (^He,^ He) interaction. 



At present, the reaction rates are generally considered to be reasonably well under control, and in the 
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CES AM code in particular, we use th e most recent estimates for the ^Be(p, y)^B (lJunghans et al. and 
''*N(p, 7)'^0 (IFormicola et al.ll2004l) . That last cross section leads to a consequent reduction of the CNO 
contribution to the nuclear energy in main sequence. Some uncertainties could remain on the screening 
effect of CNO processes but they play a rather small role in the present confrontation of the models to 
the two mentioned probes. Tables \T\ and |2] summarize our previous and present works. For a comparison 
between the values o f table 1 and the work of Bahcall and collaborators along years, the reader is refered to 



Turck-Chiezd (|2004[) . Table 2 takes i nto account the neutrino oscilla tions determined by Arn^j 



, 7 IQ-^eV^ 

and tg^f?i2 = 0.45 recommended by iBahcall & Pena-Garayl (12004) . these values are compatible with the 
most recent estimates within the error bars. 



2.1.2. The impact of the solar composition on the radiative zone properties 



One important question has emerged from the new spectroscopic estimates of the carbon, nitrogen and 
oxygen atmospheric abundances, leading to a reduction by 20-30 % of the related abundances: " Are we so 
sure that we understand properly the inner composition of the Sun in the radiative zone" ? This question has 
an important issue because the solar abundance is choosen as a reference for the stellar evolution models. 
For instance, up to recently, the abundance of these nuclides was larger than that of any other heavy element 
by about a factor 10 in Population II star s, so the change in the solar CNO abundance has a strong impact on 
the evaluation of their chemical pattern (ITurck-Chieze et al.ll2004ar). The situation is now clarified because 



the three independent analyses (IHolweger 



2001 



Asplund et al.ll2009l : ICaffau et al.ll2008b converge to about 



the same reduction. A reduction by a similar factor (30%) of th e iron content also appeared in the early 
nineties, but with rather small effects on the whole radiative zone (|Turck-Chieze & Lopeslll993h . 



As fai- as CNO elements ai^e concerned, this substantial reduction has a real impact on the radiative 
sound speed profile mainly due to the subsequent change in the oxygen opacity (see Fig. [T|) and modifies 



Bahcall et al. 



2005). The confrontation 



the depth of the convective envelope (iTurck-Chieze et al.ll2004al : 
between helioseismic indicators and results of standard solar models was examined extensively (Turck- 
Chieze et al. 2004; Bahcall et al. 2005, Guzik et al. 2005). It is frustrating to note that the sound speed 
and density profiles of the model including the updated solar chemical composition are very similar to those 
corresponding to the case where we omitted to account for the slow gravitational settling of the nuclides 



(Thouletal 



19941 : iBrun et all Il998h and that this last update destroys the apparent agreement with the 



observed one. We may thus quantify the effect of the change of composition since the inclusion of the 
microscopic processes corresponds to a reduction by 10-12% of all the heavy elements. Consequently a 
20% variation of CNO abundances is almost equivalent to the 7-9 % variation of all the heavy elements 
due to gravitati onal settling, or to an eq uivalent change of opacity coefficients if all the elements were 
concerned. See lTurck-Chieze et al.l (119971) where the impact of a lot of changes are considered and Table [T] 
for the consequence on neutrino predictions. 



^The astrophysical S- factor determined by Junghans et al.l ( 120051) is Si7(0) 
value is still discussed by Igamov & Yamiukhamedo\T ^ 20081) 



21.4 ± 0.5 (stat) ± 0.6 (syst). Some slightly higher 



-7- 



EfFectively, a change in the CNO abundances also affects the mean total opacity of the stellar plasma, 
and thus the evolution of the Sun up to its present status as well as its present structure. The knowledge of the 
opacity coefficients used in a solar model is up to now purely theoretical. The conditions of temperature and 
density in the radiative zone ensure that the plasma is only fully ionized for its main constituents, hydrogen 
and helium, but heavier species such as iron down to oxygen become partially ionized in the radiative zone 
as summarized in Fig. [U It is why the bound-bound and the bound-free interaction of photons with matter 
limit the ev acuation of the energy p roduced within the inner 25% of the radius, corresponding to half the 
solar mass (llglesias & Rogerslll996h . The small amount of iron (only 2.8 10~^ of the hydrogen contribution 
in f raction number) contributes to about one fifth of the opacity cross section in central conditions (see Fig. \T\ 



and iTurck-Chieze et al.l (119931) for more details). The oxygen is the second contributor and plays a crucial 
role to determine the base of the convective zone. 

Up to now the details of the ion interactions have been verified only indirectly through acoustic pul- 
sation eigenmodes. These modes probe plasma properties throughout the Sun from 6% Rg up to the limit 
of the convective zone. However they depend on both the detailed composition and opacity. To disentangle 
both contributions would require a high radial resolution of the seismic data allowing to capt ure some spe- 
cific b ound-bound effect through bumps in the radial sound speed profile in the radiative zone (ITurck-Chieze 
1998h . that is still out of reach at present. 



Some proposed solutions t o reconcile the standard model with revised CNO abundances and helio 



seismology (|Guzik et al. II2005I : iLin et al.ll2007l: 



Basu & Antia 



2008b have not been retained. Some others 



seem to favor the old composition (|Basu et al.ll2007l : ISerenelli et al.ll2009r) . We explore, here, the idea that 
the standard solar model (SSM) framework is questionable. Considering that any improvement has always 
been considered like a new informative fact and that a better determination of CNO was waited since the 



nineties (iTurck-Chieze et al.ll 19931) . the present facts suggest thr^ee reasonable solutions which could all exist 
simultaneously: 



The gravitational settling is not well constrained tod ay for CNO and heavy elements. The photo- 
spheric heliu m content deduced fro m helioseismology (IVorontsov et al. I ll99lb of Y = 0.25 ±0.03, and 
confirmed by iBasu & Antial (Il995h for OR\L equatio n of state has dem onstrated the need to take into 
account this slow atomic diffusion introduced first by lCox et all (119891) . This process leads today to a 



reduc t ion of about 10-15% of the He mas s fraction at the solar surface dChristensen-Dalsgaard et al. 
I993I : lThouletaDll994l : 



Brun et al. 



19981 : iMichaudI 120041) but we cannot check yet if the order of 



magnitude is correct also for CNO and iron. 

The opacity coefficients used in the models have never been verified by laboratory measurements 
under the physical conditions of the stellar interiors and may be a source of uncertainty. The high 
energy lasers developed in France (LIL, LMJ) and in the United States (NIF) will offer conditions 
corresponding to the solar radiative zone and might help constraining these quantities. We are prepar- 
ing experiments to measure the absorption energy spectrum of isolate d elements r mixtures in plasma 
conditions which are more and m ore useful for the stellar community (|B alley et al. II2OO7I biLoisel et al.l 



20091 : iTurck-Chieze et al.ll2009r) . Laboratory opacity and equation of state experiments might benefit 
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in the coming decade from an equivalent effort to the one dedicated to nuclear reaction rate mea- 
surements. A detailed verification of the interaction between photons and plasma will contribute to 
disentangle the different processes in the deep interior of stars. 

The radiative zone of the Sun is not well reproduced by the hypotheses of the SSM, and we need 
to improve the models by the addition of dynamical processes. It is one of the reasons to develop 
the Dynamical Solar Model (DSM) which must introduce the dynamics of the tachocline, the gravity 
waves generated at the base of the convective envelope and propagating in the radiative interior, the 
presence of a potential m agnetic field and the different transport processes induced by the rotation. 
Guzik & Mussak I (120101) examine the idea of r nass loss and accretion, a nd in a coming work, we 
question the hypothesis of the energetic balance (ITurck-Chieze et al.ll2010f) . 



Such a complete model of the Sun is not ready to be compared to seismic and neutrino probes, it is 
why we have developed an intermediate model, the seismic model described in the following section. A first 
step toward the DSM is then to confront all the probes with solar models including rotation, which is the 
purpose of the present paper. We shall produce models in the following sections which can be considered as 
the DSMl step. 



2.2. Prediction of neutrino fluxes, gravity mode frequencies and structural profiles from SSeM 



The quality of the seismic observations is such that one can build a seismic model. This model is 
designed to reproduce the measured solar sound speed using the classical stellar evolution equations and 
slightly changing some recommended specific absolute values of physical inputs within their own uncer- 
tainties (|Turck-Chieze et al.ll2001al . l2004al : ICouvidat et al.ll2003ah . It can be established for the new solar 
composition and will converge to the same predictions by definition. 

The interest of such a model comes from the idea that the standard framework is obviously too crude to 
reproduce all the existing observed quantities. This model is useful to improve the prediction of the neutrino 
fluxes (Tables [H and 111) or the gravity rnode frequencies because it implicitly inc ludes the known sound 
speed profile (|Turck-Chieze et al.ll2004al : Icouvidat et allboOSal : iMathur eTaPboO?! '). 



Table [T] illustrates the time evolution of the predicted neutrino flux which depends strongly on the 
central temperature and consequently, on the details of the plasma properties. The seismic model prediction 
(??) agrees remarkably well with the measured value obtained with the SNO detector (filled with heavy 
water), which is sensitive to all the neutrino flavours. For the gallium or water detector predictions, one 
needs to inject the energy dependent reduction factor due to the fact that the electronic neutrinos ai^e partially 
transformed into muon or tau neutrino s. This property of neutrinos has been confirmed last yeai^ by the 
Borexino results (|Arpesella et al.l 120081) . Doi ng so, the agreement bet ween the predictions of the seismic 
model and all the detectors is extremely good (|Turck-Chieze et al I2OO5I) . Table 2 confirms that it is also true 
for the beryllium neutrinos. 



Another interesting property of the SSeM is the fundamental periodicity Pq which characterizes the 
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Table 1: Time evolution of the boron neutrino flux prediction (expressed in 10^ /cm^/s) associated with the 
reaction ^Be (p, y) — > ^Be * + Vg — > ^He, for calibrated solai^ structural models. Also mentioned 
are the con^esponding central temperature T^, initial hehum content 7o and the origin of the improvements 
introduced in the corresponding solar model. 



Year 


Boron neutrino flux 


Tr 

C 


U 


Problem solved 


1988 Turck-Chieze et al. 


3.8 ± 1.1 


15.6 


0.276 


SSM CNO opacity, ^Be{p,y) 


1993 Turck-Chieze & Lopes 


4.4 ± 1.1 


15.43 


0.271 


SSM Fe opacity, screening 


1998 Brun et al. 


4.82 


15.67 


0.273 


SSM Microscopic difi'usion 


1999 Brun et al. 


4.82 


15.68 


0.272 


SSM Turbulence in tachocline 


2001 Turck-Chieze et al. 


4.98 ± 0.73 


15.74 


0.276 


SSeM 


2003 Couvidat et al. 


5.07 ± 0.76 


15.75 


0.277 


SSeM +magnetic field 


2004 Turck-Chieze et al. 


3.91 ± 1.1 


15.54 


0.262 


SSM Asplund composition 


2004 Turck-Chieze et al. 


3.98 ± 1.1 


15.54 


0.262 


SSM Holweger composition 


2004 Turck-Chieze et al. 


5.31 ±0.6 


15.75 


0.277 


SSeM & new 'Be(p,y), ^'*N(p,y) 


Present work 


5.31 ±0.6 


15.75 


0.277 


SSeM 




4.21 


15.51 


0.262 


SSM Asplund 2009 




5.09 


15.64 


0.273 


SSM GN and DSMl GN with slow rotation 




4.52 


15.54 


0.269 


DSMl GN moderate rotation 



Table 2: Comparison of the calibrated Solar Seismic Model predictions and the detection of the difi'erent 
neutrinos. 

SNO predictions 5.31 ± 0.6 lO'' cm"^ s'^ 
SNO results 5.21 ± 0.27 (stat) ± 0.38 (syst) lO^cm'^s'^ 
Gallium prediction without neutrino oscillations: 123.4 ± 8.2 SNU 
Gallium prediction with neutrino oscillations: 66.65 ± 4.4 SNU 
Galhum detection: GNO-^Gallex: 70.8± 4.5± 3.8 SNU; SAGE: 70.8± 5.3± 3.7 SNU 
Chlorine prediction without neutrino oscillations 7.6 ±1.1 SNU 
Chlorine prediction with neutrino oscillations: 2.56 ± 0.23 SNU 

Chlorine detection: 2.56 ± 0.16(stat) ± 0.16(syst) SNU 
^Be prediction without neutrino oscillations: 4.72 lO^cm^^s"^ 
^Be prediction with neutrino oscillations: 3.045 lO^cm^^s"^ 
BOREXINO detection 49 ± 3 (stat) ± 4 (sys) cts/day/100 tons or 3.36 ± 0.365 lO^cm'^s"! 
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Fig. 1. — Relative contribution of the most important heavy element contributors to the total opacity when one 
considers the solar composition proposed by Grevesse & Noels (1993) (left) and Asplund et al. (2004) (right) . An 
intermediate value of the oxygen contribution is now more compatible with the recent studies ( Asplund et al. 20091 



Caffau et al 



2008 



20091). We have noticed that slight differences exist between the different opacity tables due to 



slight different states of ionization of the different species (ITurck-Chieze et al 



Table 3: Characteristics of the rotating solar models computed with CESAM (Ac, B^) and STAREVOL 
(As, A's, 65,05) codes. All STAREVOL models where computed using aMLT = 1.7378 and initial helium 
mass fraction = 0.280. In this study, all the CESAM models are calibrated in luminosity and radius at 
the present age. 



Model 


D/, 


^initial 






Ac 


Mathis et al. (2004) 


3.27 10"^^ 


2.19km.s * 


2.13 km.s-^ 


A5 


Mathis et al. (2004) 


4.84 10"^^ 


2.15 km.s"' 


2.05km.s-i 


As 


Zahn(1992) 


4.84 10^^ 


2.15 km.s"* 


2.03km.s-' 


Be 


Mathis et al. (2004) 


1.10 10'^'' 


19.6 km.s"' 


2.09 km.s"! 


B5 


Mathis et al. (2004) 


3.88 10'*'' 


19.7 km.s-i 


3.08km.s-' 


C5 


Mathis et al. (2004) 


8.74 10^'' 


53.2 km.s'' 


3.03km.s-* 
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asymptotic behaviour of the gravity modes nearly equally spaced in period for frequencies below ~ 100 //Hz. 
Pq = 2n^ (^j^' 7''^^) where A'^ is the Brunt- Vaisala frequency, is the radius at the convective bottom. Be- 
fore the launch o f SoHO, there was a great dispersion on the theoretical predicted values of Pq. Following 
Hilletal.l(ll99lk its value was varying from 29 minutes to 63 minutes depending on the models. Today, 



standard and s eismic models agre e within 1 minute, which helps to predict the general properties of the 



gravity modes (|Mathur et al.ll2007h because models have been improved by the understanding of the prop- 
erty of the modes and their sensitivity. Nevertheless the present standard model, including the new CNO 
composition, largely deviates from the seismic observations, consequently the value of Pq might be better 
determined by the seismic model or any coming model in agreement with seismic observations. Of course 
the detection of individual gravity modes might improved the density profile in the core and consequently 
will put more pressure on the seismic model in the core. 

Despite the remarkable agreement between the two probes of the central region of the Sun (neutrinos 
and helioseismology), SSeM is not a physical model and SSM predictions with the new solar composition 
encounter a series of difficulties. In the following section, we present a first step to transform SSM or SSeM 
into a more realistic DSMl, by introducing rotation and the associated transport of matter and angular 
momentum in the computations. 



HOW TO MODEL THE INTERNAL ROTATION IN ID STELLAR EVOLUTION CODES 
3.1. Meridional circulation and shear-induced turbulence : adopted formalism 



The present Sun is a slowly rotating star (- 2.2 km/s at the equator, 1.7 km/s neai^ the poles). It is 
thus a reasonable first order approximation to treat rotation with a perturbative approach and to adopt the 
hypothesis of spherical symmetry in the standard model framework. We may nevertheless recall that the 
solar superficial defor mations can be measured and that the quadrupolar mo ment is of: the or der of 1.84 10"^ 
Lvdon & Sofialll995l) and will be improved by the microsatellite PICARD dThuillieilboOsI) . 



Even if departures from spherical symmetry can be considered small in the solai^ case, they ai^e suffi- 
cient for rotation to induce large-scale circulations in the ra diative and convective in terior, that will simul- 
taneously advect angular momentum, and chemical species (Busse 1982 : Zahnlll992|) . Moreover, in case of 
radial differential rotation, which could be the case in the solar deepest interior, different hydrodynamical 
instabilities may develop that will generate hydrodynamical turbulence in the radiative regions. La boratory 
Couette-Taylor experiments conducted for Reynolds numbers Re ^ 10^ by iRichard & ZahnI (119991) indicate 
the development of turbulence in the flows. This gives a good hint that in stellar interiors, where the typical 
Reynolds number is of about 10^^ (solar case), the hydrodynamical instabilities associated with differential 
rotation eventually become turbulent. 



Following IZahnI (|l992h : iMaeder & ZahnI (|l998h : iMathis & ZahnI (|2004l) . the transport of angular mo- 
mentum and chemical species by meridional circulation and shear-induced turbulence has been added to 
the stellar structure equations in order to estimate the effects of such a transport on the evolution and the 
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structure of the Sun up to the present time. 

Under the assumption of shellular rotation ensured by a strong anisotropy of the turbulent diffusivities 
D/, >> Dy, all the relevant variables and vectorial fields m ay be described as the sum of th e mean value 
over an isobar and of a second (or fourth) order perturbation ( Zahn|[l992 : Mathis & Zahnll2004[) . This means 
that the variables are projected on a basis of Legendre polynomials, an approach very similar to that used 
in helioseismology. This allows to separate angular and radial parts and thus to account for the rotational 
transport of angular momentum and chemicals over secular timescales in 1-D stellar evolution codes, while 
2-D or 3-D secular stellar evolution is yet far from maturity. 

Using this formalism, the momentum equation 



dV / ^ ^\ -> 



-VP- V0 + V.||t|| 



(1) 



becomes, when averaging over the isobai^ and using an azimuthal projection : 



5r^ dr 



dr 



4^0 



(2) 



where Q = is the mean angular velocity on the isobar of radius r and Vv is the vertical turbulent viscosity 
associated to the shear instability (see Eq.[S]below). 
C/2 is the vertical meridional velocity component : 



pgC„r[Vad-v + (0M)VJ\M* 



(3) 



with P is the pressure, Cp the specific heat, and are terms depending on the Q- and [i- distribu- 
tions respectively, up to the third order derivatives and on various thermodynamic quantities (see details in 
Maeder &Zahnlll998h . 



Equation |2] is split into a system of four first-order equations in Q. complemented by an equation de- 
scribing the evolution of the horizontal mean molecular weight fluctuations due to the vertical /i gradient for 
a given horizontal turbulence and a vertical velocity: 



-T + U2{r)^ = -6DhA 
at or 

where A = ^> and Dh is the horizontal turbulent diffusivity. 



(4) 



The equation for the transport of nuc lides by meridional circu lation and shear-induced turbulence can 
be written as a diffusion equation (see e.g. lChaboyer & Zahnlll992l) : 



dXi 
"dT 



M, 



dX 

{Anr^pf (D,, + Deff) '- 



dXi 
17 



-I- 



nucl 



dXi 

IT 



(5) 
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where X, is the mass fraction of the i^'' nuclide, the second and third terms are respectively the nuclear (nucl) 
and gravitational settling (micro) terms. D,, = Vy is the vertical component of the turbulent difFusivity. Deff 
is the diffusion coefficient associated to the action of the meridional circulation, and it is given by 

D.--'4^. (6) 



30Dft 

These transport equations have to be solved at each evolutionary timestep, adding a set of five coupled 
non-linear equations to the usual set of structure equations. T his explains the smal l number of stellar evo- 
lution codes inc luding this treatment: th e STAREVOL code, (|Palacios et al.l 120061) and the Geneva stellar 
evolution code (lEggenberger et al.l 120081) . It has been recently included in the CESAM stellai^ evolution 
code. 

Following the effects of rotation according to the above described formalism is even more difficult in 
the solar case due to the small perturbation they cause on the solar structure. In order to test the quality 
of our results, we have decided to confront the results obtained with two different stellar evolution codes 
using distinct numerical approaches. We have computed rotating solar models with both STAREVOL and 
CESAM codes. In the following, we present the codes, the numerical implementation of the Maeder & 
Zahn formalism and the main characteristics of the solar models that will be discussed in the forthcoming 
sections. 



3.2. Stellar evolution codes 

3.2.1. CESAM 

The CESAM code has been developed by a consortium of french astrophysicists to get a stellar evolu- 
tion code of high accuracy for numerous seismic uses. This code is robust and used by a large international 
community. It is free of access for the basic use ("http://www.oca.eu/cesam), that rn e ans fo r standard physics 



excluding the dynamical processes. Details of the code are described in iMorell (119971) . This code does 
not use the natural variables to avoid the singularities at the center. Moreover CESAM does not solve the 
non linear problem of the limited conditions in the physical space but in the B -splines space. This allows 
a good continuity of the functions and their derivatives. In fact, this code has been especially built for an 
helio and asteroseismic perspective, so all the variables used to calculate the frequencies of the acoustic and 
gravity modes are correctly calculated. The treatment of the transport equations by rotation has been re- 
cently improved by Marques (2010), it is now done in the physical space with a projection of the quantities, 
associated to the rotation, on the vectorial spherical harmonics. These equations are solved independently 
of the structural equations, and then all ffie interpolations from one model to the next one (in time) are done 
on the B-spline basis. 

CESAM is used by our team since 1995 (see the references in Table 1), in replacement of the Paczynski 
BINARY code used in the eighty for the predictions of helio and asteroseismic probes and for neutrinos. 
With time, the code has been enriched by the different updated pp and CNO cross sections of the hydrogen 
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buming. For the Sun, we use the recommended values of the Seattle meeting (lAdelberger et al.lll998h . then 
the most rec ent updated quoted in section 2. We have included the appropri ate screening and u se the Mitler 
prescription (IDzitko et alJI 19951) and the microscopic diffusion described in iBrun et alj (119981). CESAM is 
regularly updated with the most recent equation of state and opacity coefficients dlglesias & Rogers! 1 19961 : 



Iglesias & Roselll99q) . It included the different upd ates on the solar composition with an adapted low 



temperature opacity table (|Turck-Chieze et alJl2004ah . Th e models of this paper use the Holf atmosphere 
(see the difference between Kurucz and Holf atmosphere in lCouvidat et all (|2003ah ). 

For a comparison with the STAREVOL code, we have used the Grevesse and Noels (1993) composi- 
tion. This choice has no incidence on the results because we show mainly relative comparison and because 
the solar structure is only slightly modified by the effects of rotation. 



3.2.2. STAREVOL 



The code STAREVOL has been applied to model stars all over the Hertzsprung-Russell diagram, 
and more specifically to study the effect of dynam i cal processes on stellar evoluti o n and nucleosynthe- 



sis (IPalacios et al 



Decressin et al 



2003: Charbonnel & Talon 



20051 : iDecressin & Charbonnell l2005l : iPalacios et all l2006l : 



2009 ). Let us briefly recall he re the main ingredients used for the computation of the solar 



model. We adopt the lGrevesse & Noelsl (119931) as reference for the solar chemical composition. The nuclear 
network includes 52 species up to ^^^Cl. The associated nuclear reaction rates for charged and neutron parti- 
cle captures as well as beta dec ays have bee n kept up-to-date using the NETwork GENerator tool available 



OPAL tables (|lglesias & Rogers 



at http://astropcO.ulb.ac.be (see lSiessI (120061) for details) . At Iq w-temperature {T < 8000 K), the atomic and 
molecular- op acities ai'e given by Alexander & FergusonI (|l994l) . Between 8000 < T < 5 x 10^ K, we use the 



1996f) . The neutrino energy loss rates aie computed according to 



ItohetaL 



(Il996h and take into account the effects of plasma, pair, bremsstrahlung, recombination, and photo neutrino 
emission. 

The equation of state (EOS) is bas ed on the principle of Helmholtz free energy minimization and is 



described in detail in Siess et al. 



(|2000l) . The atmosphere is treated in the grey approximation and inte- 
grated up to an optical depth r ^ 5 x 10~^. Convection is modeled following the mixing length formalism, 
and a common parameter (Imlt - 1-7378 has been adopted for the thr^ee models presented here, ensuring 
that the solar luminosity and radius are reproduced with at most 0.1% eiTor (see Table 6). We apply the 
Schwarzschild criterion for the convective instability. Convective regions are assumed to undergo instata- 
neous mixing of chemicals and to rotate as rigid bodies. For the solar model, mass loss was not taken into 
account. The Maeder & Zahn (1998) formalism is applied, and the resulting set of five non-linear diff er- 



ential equati ons is solved using a Newton-Raphson numerical scheme according to iHenyey et al.l (Il964j) as 



Montmerle & Michaudr(| 19761) for the microscopic 



described in lPalacios et al.l(|2006h . The transport of nuclides by gravitat ional settling is also a ccounted for in 
all the rotating solar models. It is introduced using the ap proximation oflPaauette et all (119861) for the micro- 
scopic diffusion coefficient and the expressions given by 
diffusion velocity. 
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3.3. Assumptions introduced in the rotating models of the Sun 



The Sun is the only star where we can perform quantitative comparison of the impact of the above 
prescriptions on the sound speed and the rotation profiles which can now be deduced from helioseismic data. 
So it is interesting to discuss in details the role of the different terms and the choice for the prescriptions. 
In the following section, we shall present two types of models: one model with an extremely slow initial 
rotation which does not justify external magnetic braking (models A or slow), and two models with greater 
initial rotation and magnetic braking at the surface (models B and C, or intermediate and strong rotation 
models). The main characteristics of these models are given in Table |3] These cases allow to discuss the 
order of magnitude of the meridional circulation, the diffusion coefficients and the gradient of the internal 
angular velocity profile, together with their different roles in the building of the present solar rotation profile. 



For these models we have adopted iMathis et all (120041) prescription for the horizontal component of 
the turbulent diffusivity D/,: 

Dh = ^|[J^(r'^{r\lV2 -aU 2\l 
The vertical component of the turbulent diffusivity D,, is from Talon & ZahnI (1997) : 

Ric 



(7) 



D, = 



[rdrO) 



(8) 



with Ric = 1/6 the critical Richai^dson number, Kj the thermal diffusivity, Nt and A'^ the chemical and 
thermal parts of the Brunt Vaisala frequency N'^ = Nj + N^. 



Referring back to iTalon & ZahnI (|l997n . let us recall here that for the shear instability to develop, not 
only the Richardson criterion must be satisfied, but the sheared flow shall be turbulent. This is controlled by 
the Reynolds criterion. In our models, we used a critical Reynolds number Rec = lOvmoi- When Vy < Rec, 
Dy is replaced by Vmoi in the transport equations for angular momentum and nuclides. 

The choice of slow and mild rotators on the Zero Ag e Main Sequenc e (hereafter ZAMS) refers to the 
two identified populations of young stars as pin-pointed by lBouvieii (|2009h : the slow rotators for which the 
angular momentum evolves like and the fast rotators for which the angular momentum varies like Q.. 

Models A explore the extreme case of slow rotators on the Pre-Main Sequence (hereafter PMS) and 
on the ZAMS, with a surface equatorial velocity of about 2.2 km.s"' remaining almost constant from the 
ZAMS to the present age of the Sun. No braking is applied at the surface. These models do not actually 
reproduce observed stars, and can be considered more as academic cases used to show the respective impact 
of the meridional circulation and shear turbulence on the shape of the angular velocity profile even with such 
small velocities. They are also used to estimate interesting quantities relevant for the present observations. 

Models B and C are mild rotators at their arrival on the ZAMS (approximately 20 km.s"' and 50 km.s"^ 
respectively). These models undergo magnetic b raking at the arrival on the main sequence, similar to the 
solar-mass models by lTalon & Charbonnell (l2005b . 

Rotation is included from the PMS in aU our models. We assume solid body rotation as the initial state in 
the completely convective PMS stai\ When magnetic braking is applied at the arrival on the main sequence. 
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we use the formalism developed by iKawalen (119881) (see also iBouvier et all 1 19971 : IPalacios et alj|2003h in 



order to obtain the solar surface equatorial velocity when the model reaches 4.6 Gyrs 

JJ\ ..^,IR\"'IM-'-"' 



-l^-Knnt,,^-! (a>a„) <,0) 

This formulation corresponds to a field geometry intermediate between a dipolar and a radial field 
(Kawaler 1988), and it is widely used in the literature (Chaboyer et al. 1995; Krishnamurti et al. 1997; 
Bouvier et al. 1997; Sills & Pinsonneault 2000). The parameter K in Kawaler's law is related to the mag- 
nitude of the magnetic field strength, and is adjusted accordingly with the adopted initial rotation velocity. 

sat expresses the fact that the magnetic field generation saturates beyond a certain evolutionary point. This 
saturation is actually required in order to retain a sufficient amount of fast rotators in young clusters, as 
originally suggested by Stauffer & Hartmann (1987). Here we adopt Q„„ = 14Q.Q following Bouvier et al. 
(1997). Table |3] presents the main characteristics of models A, B and C computed with STAREVOL and 
CESAM. 



4. THE ROTATING SOLAR MODELS and THEIR SENSITIVITY TO THE USED 

PRESCRIPTIONS 

The computation of solar rotating models is tricky due to the subtle changes induced by rotational 
mixing on the internal structure of the Sun. The aforementioned formalism being newly introduced in the 
CESAM code, we compare first in this section the results obtained with STAREVOL and CESAM in order to 
achieve a mutual validation of the results. The numerical approach to solve the stellar structure equations in 
CESAM and STAREVOL is significantly different and also is the numerical implementation of the transport 
of angular momentum equation. We thus do not expect a perfect match between the rotating solar models 
computed with each code, but orders of magnitudes and shapes for the different profiles are expected to be 
similar. 



4.1. An extreme case: the Sun, a slow rotator on the ZAMS 

We consider first the extreme case where the initial angular momentum is small so that the model 
arrives on the ZAMS as a slow rotator (i/zams ~ 2.2km.s"') without undergoing any magnetic braking. 

Table |4] shows the evolution of several quantities among them the surface velocity in the center and at 
the surface of model Ac. With an initial rotation of about 0.1 km.s"^, the surface velocity is seen to increase 
to reach 2.19 km.s"^ at the arrival on the ZAMS (t = 390 Myrs), and then remains almost constant up to 
the present age of the Sun. In the central regions, the increase of the rotation velocity is maintained during 




Fig. 2. — Comparison of angular velocity profiles between the two codes CESAM (C) and STAREVOL (S) 
for slow rotating models Ac (solid lines) and A5 (dashed lines) at 2 Gyrs and 4.6 Gyrs. 



Table 4: Characteristics of the calibrated models computed with CESAM for the slow initial rotation 



age (Gyrs) 


R(Ro) 


L(Lo) 


Eg (%) 


Epp (%) 


EcNO (%) 


Tc 


Pc 


Vc,eq (km/s) 


Us (km/s) 





17.2 


73. 


100 








4.83 10^ 


2.1 10-3 


0.08 


0.08 


0.00001 


10.5 


32.9 


100 








7.69 10^ 


8.5 10-3 




0.12 


0.0001 


4.90 


8.7 


100 








1.6 10^ 


7.8 10-2 




0.25 


0.001 


2.20 


1.8 


100 








3.53 10^ 


8.2 10-1 




0.55 


0.015 


1.12 


0.46 


100 








6.67 10^ 


1.23 10^ 


2.41 


1.18 


0.020 


1.02 


0.63 


94 


6 





9.26 10^ 


4.39 10^ 




1.64 


0.030 


0.96 


0.89 


33 


65 





1.32 10^ 


8.01 10^ 




2.06 


0.039 


0.88 


0.71 


33 


92 


8 


1.36 10^ 


7.88 10^ 




2.19 


0.050 


0.87 


0.70 





94 


6 


13.5 10^ 


7.89 10^ 


5.03 


2.18 


0.147 


0.88 


0.72 





98 


2 


13.5 10^ 


8.26 10^ 




2.18 


1 


0.9 


0.77 





100 





13.8 10^ 


91.2 


5.6 


2.177 


2 


0.923 


0.82 





100 





14.2 10^ 


102.6 
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0.948 
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100 





14.7 10^ 


111.7 
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2.16 


4.6 


0.997 
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147.4 
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Fig. 3. — Meridional circulation velocity profiles for solar models at 4.6 Gyrs computed with CES AM (solid 
lines) and STAREVOL (dashed Unes) for very slow initial rotation (left) and mild initial rotation (right). 




Fig. 4. — Diffusion coefficient profiles for solar models at 4.6 Gyrs computed with CESAM and 
STAREVOL. The horizontal turbulent diffusion coefficient D/,, the vertical turbulent diffusion coefficient 
Dy and the effective diffusion coefficient Dgff are presented for models Ac (solid) and As (dashed) on the 
left panel and for models Be (solid) and Bj (dashed) on the right panel. 
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the main sequence, so that at the present age, vq is about 4 times larger than v^. This velocity gradient 
clearly appears in Fig. |2] and on the left panel of Fig. |7J where the evolution of the angular velocity Q. is 
displayed for models Ac and As. We may note here the very good agreement between the profiles obtained 
with CESAM and STAREVOL for solar models calibrated within 1%, thus validating the numerics of both 
codes. The building of the Q-gradient is due on the PMS to the rapid contraction of the central regions when 
radiation becomes the dominant energy transport process in the core, as well as to the combined effects of 
meridional circulation and shear-induced turbulence. 

Despite the very slow rotation of these models, a small temperature gradient is generated between 
the pole and the equator and a meridional circulation develops in the radiative interior. The velocity of 
the meridional currents, U2 is of the order of 10~^cm.s~^ in absolute value at 4.6 Gyrs as can be seen in 
Fig.[3]left. This flow is extremely small compared to the meridional circulation in the convective zone. The 
circulation consists in a single counterclockwise loop, e.g. U2 negative. The flow is directed downward, 
peaking near the base of the convective zone, extracts angular momentum from the central region to the 
base of the convective envelope. The profiles obtained with both codes are of similar shape and amplitude. 
They present irregularities mainly due to the mean moleculai^ weight variations A = ^ that introduce non- 
linearities in the system because some noise is coming from the mean molecular weight gradient. Let us 
here recall that the transport of angular momentum introduced in the codes is treated in a self-consistent 
but simplified manner. Such models give hints, order of magnitudes and global shapes, but should not be 
expected to deliver exact profiles. 

The diff'usion coefficient profiles are also similar in the two codes (see Fig. |4] left). The effective 
diffusion coefficient D^yy is proportional to {1/2)^ /D^ and consequently is very small (about 10"^ to 10"^ 
cm^.s"') due to the slow meridion al circulation veloci ty and the large horizontal component of the turbulent 



diffusivity ensured by the use of iMathis et all (12004]) prescription. The Reynolds number of the sheared 



flow is lower than the critical Reynolds number Rec oc lOvmoi in the bulk of the radiative zone. The flow 
is not turbulent, so the vertical component of the turbulent diffusivity D,, is essentially set equal to the 
local molecular viscosity (see § 13.31). The eff"ectiy e diffusion coefficient is actually much smaller than the 



microscopic diffusivity (see e.g. iBrun et al.lll999h . which is of the order of 10 cm^.s"^ in solar models. 
Consequently the microscopic diffusion is the most efficient process to transport chemical species below the 
convective envelope in the models with slow rotation. This shows up in the helium mass fraction profile 
presented in quadrant III of Fig. [5] There, the step in the helium profile at the base of the convective zone 
reveals the gravitational settling of helium in this region, and this, regardless of the existence of meridional 
flows and shear. If the young Sun was a slow rotator, the chemical stratification probed by helioseismic data 
in the region below the convective envelope would be the same as if the Sun had not rotated (provided that 
no other dynamical process than rotation is included). 

Of course, one cannot ensure that the diffusivity coefficients are properly estimated, so it is interesting to see 
the impact of a dif ferent prescription for them. Using the STAREVOL code, we compute model A's using 



the prescription of lZahnI (11992!) for D/,: 



DH^r[2V2-aU2] 



(11) 
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Figure|5]shows the profiles of the angular velocity Q, the meridional circulation velocity U2, the helium 
mass fraction Y and the diffusion coefficients D/,, D,, and Dgflf for models A5 and A' 5. As expected, the 
horizontal turbulent diffusion coefficien t D/, is three orders of magnitudes smaller when using Zahn's 1992 
prescription (see also lMathis et al. Il2004l) . This translates into an increase by the same amount of the effective 
diffusivity Dgff . However, Dgff remains smaller than D,, and Vi^jc in model A'5 , so that the helium and angular 
velocity profiles reamin unchanged compared to model A5 . From this comparison, we may conclude that 
in the case of a Sun that has always been rotating slowly, the choice of the prescription for the horizontal 
component of the shear-induced turbulent diffusivity does not affect significantly the structure, rotation and 
initial chemical stratification. 



4.2. The Sun as a mild rotator on the ZAMS 

We have also computed two models with a higher initial angular momentum, corresponding to a surface 
velocity of about 20 km.s"' and 50 km.s"' at the arrival on the ZAMS. These models undergo a strong 
angular- momentum extraction on the early main sequence to reach a surface velocity of about 2 km.s"' 
at the age of the present Sun. This phenomenon associated to the idea of a magnetic braking due to the 
decoupling from the disk environment is modelled using the Kawaler formula as detailed in the previous 
section. 

Fig. |6] compares the results of the two models Bj and Cj at 4.6 Gyrs. The diffusion coefficients D,, and 
Deft entering the diffusion equation for the chemical species are similar in models and C^, so they 
lead to a same profile for the helium mass fraction. The angular- velocity is slightly larger in the central 
regions in model . This model is submitted to a more efficient braking in order to reach a surface velocity 
of about 2-2.2 km.s"^ at the age of the Sun. In fact, when magnetic braking is initiated, the meridional 
circulation velocity U2 is five times smaller in model C5 than that found in model B5 in absolute value, 
but it is positive (negative in model B^) generating an efficient inward transport of angular momentum 
and leading to a faster and larger increase of the centra l velocity than in model B ,-^. The behaviour found 



for model C.y is in all points similar to that reported by lEggenberger et all (120051) . and also resembles that 



described by iMeynet & Maedeii (120001) for more massive and fast rotating stars. The evolution further on 
the main sequence is similar for models B^ and C5, with U2 rapidly becoming negative in the later (see 
figure top right of Fig. O. The profiles displayed in Fig.[6]show no significant difference for the diffusivities, 
meridionla circulation and helium profile at 4.6 Gyrs, so that in the following we will focus on models B5 
and model Be. 

The meridional velocity profiles at 4.6 Gyrs for models B are shown on the right panel of Fig. [3] 
The velocity of the meridional currents is 3 orders of magnitude larger compared to models A5 and Ac 
previously described. At the ZAMS, the meridional circulation is maximal below the convective envelope 
as a result of the strong magnetic torque appUed at the surface. U2 is then 2 orders of magnitude larger than 
it will be later on the main sequence, when the efficient spin down of the surface layers is accompanied by a 
global weakening of the meridional circulation. The form of the meridional circulation velocity profile U2 
is maintained over the main sequence evolution, and mainly consists in one counterclockwise cell peaking 
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Fig. 5. — Profiles of the angular velocity O, the meridional circulation velocity U2, the helium mass fraction 
Y and the diffusion coefficients for the solar rotating models As (dotted lines) and A'5 (solid lines) at 4.6 
Gyr. The hatched areas indicate the convective envelope. 
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Fig. 6. — Profiles of the angular velocity O, the meridional circulation velocity U2, the helium mass fraction 
Y and the diffusion coefficients for the solar rotating models B5 (solid lines) and Cj (dotted lines) at 4.6 
Gyr. The hatched areas indicate the convective envelope. 
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below the convective envelope and extracting angular momentum from the central regions, this tirne more 



efficiently than in the slowly rotating solar models. Similar to what was obtained by lDecressin et al.l (120091) 
for their 1.5 M© model, the meridional circulation is mainly driven by the local losses of angular momentum 
due to the magnetic braking. Meridional circulation and shear turbulence do not ensure an efficient transport 
of angular momentum. One notes that U2,max is larger in model by 25% compared to model B^. 

The gradient of angular velocity is much larger in these models compared to models Ac and As , and Q. 
increases by more than one order of magnitude from the surface to the solar core. We will discuss in more 
details the implications of such a steep profile when confronted to helioseismic data. The match between 
the two codes shown in Fig. |7] is not as good as that found for slowly rotating models. This is due to the 
difference in the numerical treatment used in both codes and the way the braking is introduced. The rotation 
profile is even steeper for model computed with CESAM, in which braking occurs very efficiently on the 
ZAMS so that the Q. profile does not evolve at all in the last 3.6 Gyrs. In model Bs, angular momentum 
is extracted from the central parts as the star evolves and both the central and surface values of the angular 
velocity decrease with increasing time so that the gradient is less steep at 4.6 Gyrs than it is at 1 Gyr. 

Diffusion coefficients are displayed in the right panel of Fig. |4]at 4.6 Gyrs. They have similar shapes 
and amplitude, which implies that the chemical stratification should also be similar in models B5 and Be- 
Although U2 is larger by 3 orders of magnitude compared to models A, the horizontal component of the tur- 
bulent diffusivity is also much larger leading to D^yy that is again smaller than 1 cm.s~\ which remains much 
smaller than both the local molecular diffusivity and the microscopic diffusivity. Meridional circulation is 
thus not contributing significantly to the transport of chemical species. On the other hand, in these mild 
rotating models, the Reynolds number associated with the sheared flow is larger than the critical Reynolds 
number Rec, and the flow eventually becomes turbulent. The vertical turbulent diffusivity reaches values 
of the order of 100-1000 cm.s"^ below the convective envelope, which is at least 10 times larger than the 
microscopic diffusivity in this region. As a consequence, the abundance profiles are flattened in both models 
B5 and Be (see e.g. see solid line in quadrant III of Fig. |6]for model B5), and differ thus significantly 
from those derived for standard solar models. If the young Sun was a fast rotator and experienced magnetic 
braking during the early main sequence without any inhibition of the action of the shear-induced turbulent 
mixing, the chemical stratification of the present Sun below the convective zone might bear the signature of 
this past rotational history. 



5. CONFRONTATION BETWEEN THEORETICAL AND OBSERVATIONAL PROFILES 

5.1. The solar observed rotation profile 

After more than 10 years of observations with GOLF and MDI instruments located onboard SoHO, 
and with the facilities of the GONG groundbased network, we have derived very important constraints on 
the radiative zone rotation profile. The determination of the splittings of a large number of acoustic modes 
has definitely established that the rotation profile in the part of the radiative zone which is not influenced 
by the nuclear reaction rates, e.g. the region between 0.25 to 0.68 Rq, is really flat with invisible latitudinal 
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Fig. 7. — Evolution of the angular velocity profiles at six diff"erent ages as indicated on the plots, for models 
Ac and A5 (left panels) and Be and Bj (right panels). The age of 45 Myr con^esponds to the arrival on the 
main sequence. 
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difFerential variation dCouvidat et alj|2003bl : EfF-Darwich et al.ll2008r) . in great contrast with the latitudinal 
differential profile of the convective zone. In the central region, below 0.25 Rq, only gravity modes may 
inform on the rotation profile. Some individual candidate modes at high frequency have been observed by 
the GOLF instrument dTurck-Chieze et al.ll2004bl : lMathur et al.ll2007l : lGarcia et allboosb . and two solutions 
have been extracted from these first data depending on the interpretation of the pattern attributed to an ^ = 
2, n - 2 mode: a slightly reduced rotation rate or an increase by about a factor 3 in the central region 
(|Mathur et aLlbOOSh . Nevertheless, only the solution of a rapid rotating core is compatible with the other 
kind o f detection using the asym ptotic behaviour at low frequency which is detected with a very high proba- 
bility ( Garcia et al.ll2007 . 2008 ). It could presume a rather complex solar core rotation larger than the rest of 
the radiative zone, with a different axis, and may be with some manifestation of the radiative zone magnetic 
field. The other solution cannot be totally excluded but is not the most probable. Again, the core profile 
needs to be confirmed with extensive data from SoHO and with an improved instrument like the GOLF-NG 
concept ( Turck-Chieze et aDl2006 , 2008 ) which will measure velocity displacements at 8 heights in the so- 
lar atmosphere. Considering that the complete solar rotation profile might be accessible to the observation, 
it is important to predict it properly down to the central region. 



5.2. The theoretical solar rotation profile 

This paper is a new step toward a detailed understanding on how the different dynamical processes can 
influence the present rotation of the Sun. In this study we have followed three different cases because we 
can only have some indirect information on the solar internal rotation profile at the end of its contraction 
phase. We have already seen in the previous section that the initial angular momentum content and the 
rotation history, in particular the fact of undergoing magnetic braking on the ZAMS, shape the angular 
velocity profile, determining the absolute value as well as the gradient. In this section, we wish to focus on 
the evolution of the angular velocity profile. We present hereafter a detailed analysis of the construction of 
the predicted internal profile at the age of the Sun and confront these predictions to helioseismology. 



Compared to previous studies by lEggenberger et al.l (120051) and lTalon & Charbonnell (120051) where ro 



tation was only included from the ZAMS, with an initial flat profile, we have decided to follow rotation from 
the PMS when the star is actually completely convective. We also assume solid-body rotation as the initial 
state, but by the time the models reach the ZAMS, the convective region has retreated to the surface and the 
O-profile is not flat anymore. 

Tables |4] and [5]list the evolution of the radius, the luminosity, the central and the superficial velocity for 
models Ac and Be and some other indicators of the contraction phase and of the involved nuclear reactions. 
The evolution of the angular velocity profile for these models is also presented in the left panels of Fig. |7l 

For model Ac, where the initial rotation rate is of about 0.08 km.s^' , the contraction of the star during 
the PMS rapidly generates a radial gradient in the Q profile. The very slow meridional currents and the 
small amount of shear generated in such slowly rotating models are not very efficient to transport the angular 
momentum. The differential rotation established during the early evolution by the contraction of the inner 
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regions is maintained and slightly amplified during the main sequence by the advection term. The predicted 
overall cont rast between core and sii rface velocities is of about a factor of 4, which is comparable to that 
obtained by lEggenberger et alj (120051) in their slowly rotating model. 



Models B and C are s imilar to those p r esente d in the three main previous s tudies of the effect of r otatio n 
on the solar evolution by IChaboyer et all (|l995h : lEggenberger et al.l (120051); ICharbonnel & TalonI (120051) . 
The evolution of the angular velocity profile is in agreement with that obtained in the later study, where 
STARE VOL models were also used . C omparing the evolution di splayed in Fig. |7] with that of Fig. 1 and 
Fig. 2 of ICharbonnel & TalonI (120051) and lEggenberger et all (120051) respectively, one can measure the impact 
of the numerical approach. In all these cases, if the general forms and amplitudes of the angular velocity 
profile at the age of the Sun show encouraging similarities (see also Fig. [Hi, the in between evolution can be 
quite different. Models B5 and Be undergo a very strong braking at the arrival on the ZAMS leading to a 
quick spin down of the convective envelope and a contrast between the surface and central angular velocities 
of about a factor 15 for model B^ already at 1 Gyr, and of about 7 in model B5 at the same age. In the later 
model, the contrast increases during the main sequence to reach a factor of 10 between the core and the 
surface at 4.6 Gyrs. In both models, the Q-gradient steepens in the bulk of the radiati ve zone as the star 



evolve s. The profiles are flater in model Bj , yet not as flat as those presented in Fig. 2 of lEggenberger et al. 
J2005h . 



A detailed analysis of the relative importance of meridional circulati on versus shear-turbulence in our 
models A5, B5 and Cs using the tools presented in lPecressin et al.l (|2009r ). demonstrates that the behaviour 
of these solar models is in all points similar to that obtained for the 1.5 Mq model shown in that paper. 
The angular- momentum transport is ensured by meridional currents. In the case of models B and C, these 
currents, although slow, are primarily generated by the action of the applied torques resulting from the action 
of magnetic braking. In the case of the slow model A, no torque is applied and the angular momentum flux 
is smaller yet also dominated by the meridional circulation. The flux of angular momentum due to shear 
turbulence is more than 5 orders of magnitude smaller than that attributed to the meridional circulation in 
all the 3 cases. 



5.3. Discussion 



These detailed calculations are self-consistent and rely on the use of a set of parameters that have 
been extensively tested thro ughout the HR diagram They consider initial rot ation generally smaller than 
used in previous litterature ( Chaboyer et al. 1995 : Krishnamurthi et al. 19971) and show slower differen- 
tial rotations in the core th an found in previous works dPinsonneault et al. 1 1 19891) or the most recent of 
Charbonnel & TalonI (l2005h . Figure |8] shows a first comparison between the different profiles and the pre- 
sumed observed rotation profile. It is clear that the models with an ini tial reasonable rota tion velocity -1- 
magnetic braking are not supported by the observations. As indicated by iGarcia et al.l (|2007h . the helioseis- 
mic data reject an increase by a factor 10 for the core rotation. On the contrary the slow initial rotation leads 
to a profile for the present Sun for which the contrast between the core and the surface angular velocity 
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Table 5: Characteristics of the calibrated models computed with CESAM for an initial rotation of 20 km/s 



age (Gyrs) 


R(Ro) 


L(Lo) 


Eg (%) 


E/,p (%) 


Ecivo (%) 


Tc 


Pc 


Uceq (km/s) 


i/s (km/s) 





17. 36 


72.4 


100 








0.48 10^ 


2.05 10-3 


0.47 


0.47 


0.00001 


10.62 


32.5 


100 








7.58 10^ 


8.26 10"3 




0.74 


0.0001 


4.95 


8.6 


100 








1.58 10^ 


7.55 10-2 




1.51 


0.001 


2.28 


1.92 


100 








3.38 10*^ 


0.73 




3.18 


0.005 


1.36 


0. 63 


100 








5.32 10^ 


4.1 




5.51 


0.011 


1.12 


0. 46 


100 








6.62 10^ 


12.15 




7.99 


0.026 


1.024 


0.897 


69 


31 





11.9 10*^ 


75.53 




16.78 


0.039 


0.88 


0.71 





93 


7 


13.58 10^ 


79.18 




20.25 


0.045 


0.87 


0.69 





93 


7 


13.56 10^ 


78.68 


26.7 


19.6 


0.147 


0.88 


0.72 





98 


2 


13.5 10*^ 


82.6 


30.8 


5. 


1 


0.9 


0.76 





100 





13.8 10*^ 


91. 




2.19 


2 


0.921 


0.82 





100 





14.2 10*^ 


102.4 


38.8 


2.08 


3 


0.945 


0.88 





100 





14.7 10*^ 


116.4 


37.3 


2.086 


4.6 


0.992 


1.000 





99 


1 


15.61 10^ 


146.6 


39. 


2.09 



Table 6: Comparison of the calibrated solar models computed with the CESAM code and those of the 
STAREVOL code. X; and Y, are respectively the initial hydrogen and helium mass fractions. Xc, Y^, 
and Pc are respectively the present central hydrogen, helium, temperature and density. Y^ is the present 
superficial helium mass fraction, a is the mixing length parameter, vbzc the position of the base of the 
convective zone, (Z/X) , the mass fraction of the heavy element in considering the GN composition. 
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CESAM 


CESAM 


CESAM 


STAREVOL 


STAREVOL 


STAREVOL 


Model 


Seismic 


SSM 


slow (A) 


moder. (B) 


slow (A) 


moder. (B) 


fast (C) 


X,- 


0.7064 


0.7075 


0. 7078 


0.7117 


0.7015 


0.7015 


0.7015 


Y,- 


0.2722 


0.2729 


0.2727 


0.2690 


0.280 


0.280 
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0.3371 


0.361 
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0.3522 


0.3397 


0.3422 


0.3456 




0.6428 
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15.67 
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147.8 


147.1 


146. 
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Y, 


0.251 
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0.253 
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0.0245 


0.0245 
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Fig. 8. — Comparison between the solar internal profile predicted by the different models and that deduced 
from helioseismology. For models Bj and Cs, the actually ploted value is Q/2;r-0.2325 rather than Q./27T 
in order for the surface value they reach to match that reached by model Be- The data points down to r/R© 
= 0.2 are deduced from the acoustic mode splittings determined by the observations of GOLF, MDI and 
GONG instruments (from Eff-Darwich, 2008). The data points in the core correspond to the supposed core 
rotation extracted from the potentially observed gravity modes. These values have still large enor bai^s and 
need to be confirmed (inspired by Turck-Chieze et al., 2004b, Garcia et al. 2007 and Mathur et al. 2008). 
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is very similar to the presumed observed one. The contrast is similar but of course none of the different 
models investigated here create a flat rotation between 0.25 and 0.7 solar radius. Within the framework 
of the present study, models As and Ac corresponding to an extremely slow rotating young Sun are more 
compatible with the central observations. However, this cannot be taken at face value be cause thi s mode l 
seems rather irrealistic because it did not contain any magnetic braking (see Figure 2 of iBouvien (120091) . 
So one can imagine a slightly greater initial rotation, typically 5 or 10 km/s at the ZAMS, with a central 
rotation value slightly eroded by some other process. The discrepancy between helioseismology and the 
profiles presented here calls for the inclusion of additional dynamical processes that efficiently extract angu- 
lar momentum from the radiative interior . This result confirms the most recent works. Probably the radiative 



magnetic field (lEggenberger et al 



2005 



Yang & Bill2006l: iDenissenkov & Pinsonneaultil2007h or the inter- 



nal gravity waves dCharbonnel & Talodl2005h or both, could flatten the profile but their role could be smaller 
than thought previously. 



5.4. The sound speed profile, composition and neutrino predictions 

Figure |9] compares the squared sound speed profile difference be tween seismic observa tions and the 



models computed with CESAM. The used acoustic modes ai^e given in lCouvidat et all (12003 al) . We present 
the three models (standard, low rotation and moderate rotation) computed with the version of CESAM used 
in the present study. These models slightly differ from the standard models already published due to the 
large reduction of the (^^N, p) reaction rate mentioned in section 2. The CNO combustion is rather small. 



but this strong decrease, previously studied in iTurck-Chieze et al.l (l2001br) . leads to a negative sound speed 



difference in the core partly compensated by a small increase in the rest of the radiative zone. These two 
differences increase the difference with observations. But what we would like to emphasize in this figure is 
the intercomparison between models with and without rotation. 

Even if these models must be more representative of the real Sun, the sound speed profile obtained for 
the two calibrated CESAM rotating models A and B are in less agreement with the observed profile than 
the standard model. In fact the rotation-induced mixing has reduced the effect of the microscopic diffusion 
(see also § 4) partly inhibited by the macroscopic motions. As discussed earlier, such effect is practically 
not visible in models A, where microscopic diffusivity remains larger than both the vertical turbulent and 
the effective diffusivity, but it is clearly visible in models B and C rotating faster on the ZAMS. 

As a natural consequence, Table|6]shows also that the central temperature is slightly reduced in rotating 
models in comparison with the SSM. For the CESAM models calibrated in radius and luminosity at 10"^ 
level, the predictions for the neutrinos are always worse than those of the seismic model (see details in Table 
1). Moreover in CESAM models where we have kept superficial Z/X constant, one may notice that the 
surface helium content also increases with the initial rotation due to turbulent flow just below the convective 
zone, as previously mentioned in Section 4 for the STAREVOL results but by a smaller amount (comparison 
of the two codes in Table 6). This will contribute to reconcile the observed helium content to the predicted 
one when one will use the most recent composition of Asplund et al (2009), if the initial rotation is sufficient 



-30- 



to allow turbulent flow at the base of the convective zone and in the tachocline. 



6. CONCLUSION AND PERSPECTIVES 

In this paper we have shown the following facts: 

- We have examined three initial rotation rates (models A, B and C), choosing initial angular- momentum 
contents corresponding to 2.5 km/s, 20 km/s and 50 km/s at the ZAMS. The last two values have also been 



adopted in other studies of the rotating Sun and s olai'-type stars (lYang & Bill2006l : Eggenberger et al.ll2005l : 



Talon & Charbonnel 2005; Chaboyer et al.|[l995l) . At 4.6 Gyrs, the thr^ee cases show a radial gradient of 



rotation in the core. Its amplitude depends strongly on the initial rotation rate. If this one is small, the 
radial gradient is estabhshed during the contraction phase and is amplified during the subsequent evolution 
up to the present Sun. For models rotating faster, one notes much steeper gradients at the age of the Sun, 
the angular momentum losses associated with magnetic braking on the ZAMS are responsible for the rather 
small decrease on the ZAMS. 

- The transport of angular momentum in the solar radiative zone during the main sequence appears 
extremely small and the meridional circulation in the radiative zone is smaller by 10 orders of magnitude 
in comparison with the observed convective meridional circulation velocity at 99 % R©. This process is 
even slower than the microscopic diffusion (gravitational settling) that we use in the radiative zone. As a 
first consequence such implementation in a stellar evolution code is delicate mainly because one needs to 
solve four coupled equations with derivation of quantities that practically do not vary. In order to settle our 
conclusions concerning the form and order of magnitude of the different quantities associated to the rotation- 
induced mixing, we have confronted models obtained with two different codes, CESAM and STAREVOL, 
using distinct numerical approaches and the old composition of Grevesse and Noels. We have shown, for 
models A and B, that they lead to rather similar results and the same kind of profile for the present Sun, 
thus validating the results and numerical approaches. Such a very large difference between the meridional 
velocity in the radiative zone (10~^ - 10~^ cm/s) and in the convective zone ( m/s) would naturally produce 
some turbulent hydrodynamical layer called and this is an interesting result of the present calculation. 

- Although the combined effect of meridional circulation and shear-induced turbulent associated to 
rotation is small, this study allows us to present radial rotation profiles that can be directly compared to 
the seismically observed one. This study sustains the idea that the Sun was not a rapid rotator after the 
contraction phase. The angular velocity profile we get for models A is not far from the presumed solar 
one in the core. Of course this model should not be considered at face value since it does not take into 
account any magnetic braking generally observed in young clusters, but it can be useful to estimate the 
interplay between processes. The second model (moderate rotation) is more realistic but shows a greater 
steep gradient in disagreement with the published detection of gravity modes. So one can imagine that the 
Sun has been in an intermediate case arriving at 5-lOkm/s at ZAMS and that its rotation profile would have 
being eroded by some other process. 



- Let us stress however that other dynamical processes known to generate efficient transport of angular 
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momentum such as magnetic fields and internal gravity waves were not yet included in our models. The 
inclusion of these processes is in the scope of a further study, and is expected to significantly aff"ect the 
choice of the preferred mode l. For instance, in ternal fossil magnetic field produces certainly a very small 
effect on the solar structure (IDuez et al.ll2009al) but may lead to efficient transport of angular momentum 
that would help flattening the angular velocity profile in the radiative zone. But before we would like to 
study t he activity of the very young Sun. Previous works show how magnetic field can modify the rotation 
profile (|Eggenberger et al iDenissenkov & Pinsonneaultll2007r ) but the action could be impr oved by the 



introduction of a more sophisticated field topology which preserves the stability of such a field (IDuez et al. 



2009bh . The understanding of the magnetic field during the contraction phase has probably a crucial impact 
on the radial rotation gradient and deserves premature estimate of lithium and beryllium destruction at this 
stage. 

- Other processes may modify the present conclusions. For example in the present treatment we neglect 
the fundamental role of the tachocline which must be established at least since the arrival on the main 
sequence. The hydro (or magnetohydrodynamical) nature of such a region may alter the angular momentum 
and/or chemicals tranport due to the internal rotation but we have already noticed that a crude t reatment of 



this re gion might slightly amplify the present tendancy on chemicals gradient and structure effect (IBrun et al. 



1999h . In the present study, we note that if the initial rotation is accompanied by an efficient magnetic 
braking, it generates some turbulent flow at the base of the convective zone which smoothes the helium 
profile and increases its abundance at the surface. 

- We note that the impact of the rotation on the solar structure is rather small. As the transport of angular 
momentum and chemicals goes from the radiative zone to the convective zone, it implies a slight reduction 
of the central temperature and of the helium content in the radiative zone. Consequently it increases slightly 
the present discrepancy between model and the observed sound speed. In fact, one ca nnot exclude other 
momentum transport which may come from the convective zone and play the inverse role (IGaraud & Garaud 
20081 : lOough 1120091) . All these other processes must be included in stellar- evolution codes before getting 
a proper DSM. We see in this study that the description of the dynamics of the PMS phase (especially 
the related magnetic field evolution) will be a crucial issue too. In the comparison between CESAM and 
STARE VOL we have noticed some non negligible difference during the contraction phase for models with a 
moderate (the same for high) initial rotation rate, this phase illustrates the sensitivity of the rotation gradient 
to the numerical details and to the magnetic braking. Moreover we show in this study that independently 
of the initial rate, the central rotation value does not change by more than 50% during the main sequence. 
So the final comparison of the DSM to the observations will partly depend on the way we shall treat the 
contraction phase, the inner corresponding magnetic field evolution and the magnetic braking phase. The 
same conclusion was already reached discu ssing the problem of lithium in young stars and in the Sun 
ferun et al.ll 19991 : Ipiau & Turck-Chiezell2002l) . 



- The sound speed profile is practically unchanged when the slow initial rotation is assumed. So even 
if this model is probably more representative of the dynamics of the solar interior than the standard model, 
it is still an incomplete model and its predictive character remains limited. For this reason, and considering 
the consequent discrepancy on the sound speed predicted by the present standard model and the observed 
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one, which will be amplified by the recent CNO composition, we continue to recommend the seismic model 
for any global predictions (gravity modes or neutrinos). We have shown in this paper that the predictions of 
the seismic model are in very good agreement with all the neutrino detections including BOREXINO. 

We would like to thank PieiTC Morel for his dedicated effort to introduce the dynamical processes in 
CESAM code, Stephane Mathis and Jean Paul Zahn for very interesting discussions and the anonymous 
referee for helping us to improve the paper. 
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Fig. 9. — Radial squared sound speed and density differences between o bservations and models computed 
with the CESAM code. Seismic model ( Turck-Chieze et aPboOla , 2004a ) is in black with error bars coming 
from the seismic data. The present models (standard GN composition) without rotation, with low rotation 
(model A) or moderate rotation, of 20 km/s initial model with a magnetic break of the surface at the arrival 
on the main sequence (model B) are respectively in green (dotted dashed line), blue (large dashed Une) and 
red (small dashed line). 



